Comparative study of rainfall prediction based on different decomposition methods of VMD

Rainfall forecasting is an important means for macro-control of water resources and prevention of future disasters. In order to achieve a more accurate prediction effect, this paper analyzes the applicability of the "full decomposition" and "stepwise decomposition" of the VMD (Variational mode decomposition) algorithm to the actual prediction service; The MAVOA (Modified African Vultures Optimization Algorithm) improved by Tent chaotic mapping is selected; and the DNC (Differentiable Neural Computer), which combines the advantages of recurrent neural networks and computational processing, is applied to the forecasting. The different VMD decompositions of the MAVOA-DNC combination together with other comparative models are applied to example predictions at four sites in the Huaihe River Basin. The results show that SMFSD (Single-model Fully stepwise decomposition) is the most effective, and the average Root Mean Square Error (RMSE) of the forecasts for the four sites of SMFSD-MAVOA-DNC is 9.02, the average Mean Absolute Error (MAE) of 7.13, and the average Nash-Sutcliffe Efficiency (NSE) of 0.94. Compared with the traditional VMD full decomposition, the RMSE is reduced by 7.42, the MAE is reduced by 4.83, and the NSE is increased by 0.05; the best forecasting results are obtained compared with other coupled models.

Generally speaking, the predictability of precipitation is derived from the "memory" of atmospheric circulation 1 .However, this "memory" typically lasts only a few weeks, which is insufficient for monthly precipitation predictions.With advancements in related fields, research has discovered the existence of "long memory" in oceanic and land surface processes that interact closely with atmospheric motion, offering the potential for monthly precipitation predictions 2 .Rainfall prediction is typically categorized into statistical methods relying on correlation relationships and physical-statistical methods based on physical factors.Nevertheless, these approaches face challenges in accurately capturing the inherent nonlinearity (complexity, diversity, and instability) of water anomalies and the uncertainties affecting precipitation anomalies due to the demanding data requirements and the uncertainty in detailed regional initial condition descriptions.In contrast, data-driven machine learning methods have shown promise in precipitation prediction, thanks to their powerful nonlinear mapping capabilities 3 .These methods bypass the intricate physical boundary conditions between influencing factors and precipitation, concentrating solely on identifying the optimal mathematical relationship between them.This renders machine learning-based monthly precipitation prediction methods highly practical and convenient.Among these methods, the decomposition-reconstruction-based machine learning model for monthly precipitation sequences offers distinctive convenience as it solely necessitates historical precipitation data as input to simulate future precipitation output.Narimani et al. 4 used two methods, singular spectrum analysis (SSA) and empirical mode decomposition (EMD), to decompose the daily precipitation data, combined with light gradient-enhanced tree (LightGBM) and extreme gradient-enhanced tree (XGBoost) to make predictions.The VMD can adjust the number of modal decompositions on its own to effectively avoid modal aliasing 5 and reduce the complexity of the sequence data, which is a commonly used and efficient time-frequency analysis tool with more advanced performance in the field of extracting effective information of precipitation sequences and unfolding sequence prediction research, and it has been verified to outperform the other decomposition methods in the relevant applications 6,7 .However, most of the previous machine learning prediction models based on VMD decomposition use full decomposition to construct the samples, which may result in mixing future precipitation information in the training samples, and thus the proposed method does not have the significance of real precipitation prediction.Therefore, some "stepwise decomposition" models for real forecasting are emerging, but there is little research on this method in the field of precipitation forecasting.
In addition, the prediction performance of the model is often affected by the hyperparameters set by human beings, and the intelligent optimization algorithm is one of the main means to prefer the hyperparameters.Adaptive Particle Swarm Algorithm (APSO) proposed by Zhang et al. 8 optimizes the hyperparameters of Bidirectional Long Short-Term Memory Network (BiLSTM), obtaining higher prediction accuracy and faster convergence.Wang et al. 9 proposed a combined wavelet decomposition-prediction-reconstruction model (WDPRM), which consists of a wavelet decomposition, Particle Swarm Optimized Support Vector Machines (PSO -SVM) and BP network optimized by Artificial Bee Colony Algorithm (ABC-BP) to obtain good prediction results.Mohammed et al. 10 incorporated the Artificial Bee Colony algorithm, Particle Swarm Optimization, and Imperialist Competitive Algorithms to improve and optimize the internal parameters of the Artificial Neural Network (ANN) method for delivering a robust approach to identify non-linearities in rainfall patterns.
From the current research situation at both national and international level, the precipitation series are highly nonlinearized and noisy, and the use of data decomposition preprocessing techniques (EMD, SSA, Complementary Ensemble Empirical Mode Decomposition (CEEMD), and VMD, etc.) is a necessary means 11 .However, most of the researches directly use it for the decomposition and denoising of complete datasets (including training and test sets), and these decompositions and denoising algorithms often suffer from the "boundary effect", and the samples constructed in the training period may use the "future data information" which is unknown in reality, and the full decomposition may be difficult to satisfy the requirements of the actual forecasting; There are numerous different types of machine learning models, the traditional recurrent neural network structure is somewhat limited in the utilization of long-ago information (long memory), and the long short-term neural network (LSTM) requires a large amount of computational resources to process long-distance dependency sequences, in order to solve this problem, a variety of neural networks with an external storage mechanism have been designed, based on which Graves et al. 12 have proposed the DNC, which has been improved on the Neural Turing Machine (NTM) based on the improvement of the storage management method, and the inclusion of temporal memory links allows the DNC to jump to read or update the memory information, thus combining the advantages of recurrent neural networks and computer processing.In addition, the hyperparameters of some machine learning models greatly affect the effect of precipitation prediction, and previous studies have mostly used particle swarm algorithms 13 and genetic algorithms 14 and other algorithms with limited ability to find the optimum, and it is necessary to further investigate the related new optimization algorithms.
Hence, this study scrutinizes the shortcomings of previous full decomposition-based precipitation prediction models for the Huaihe River Basin, which are not suitable for practical forecasting.Instead, it employs novel stepwise decomposition sample construction techniques that exclude future precipitation information in assessing real monthly precipitation prediction accuracy.Additionally, it harnesses the AVOA for refining the DNC hybrid learning neural network, incorporating initial population modification and parameter combination optimization search enhancements.The study concludes by comparing the forecasting accuracy of each coupled model and analyzing the optimal decomposition forecasting method, offering valuable insights for other precipitation forecasting research.

Research methodology VMD with different decomposition methods
The rainfall time series data are relatively single, and the internal frequency information and change rule are difficult to be directly explored.Considering it as a nonlinear time series signal, the VMD algorithm is used to decompose it into a number of intrinsic mode functions (IMFs), so as to visualize the hidden information such as the inherent cyclical trend of the time series, and at the same time increase the amount of data information for the prediction model.However, in recent years, some studies in other fields have questioned the validity of this experimental framework, arguing that its true ability is "Hindcasting" rather than "Forecasting" 15,16 .The "full decomposition" approach uses the data information of the test period in the training period, which leads to the distortion of the model forecasting results.In this study, the applicability of the "full decomposition" and "stepwise decomposition" rainfall forecasting models in the Huaihe River Basin is compared.

Full decomposition of VMD
Variational mode decomposition (VMD) is a commonly used adaptive and fully recursive signal sequence processing method, which first requires the user to set two parameters-the number of decompositions K and the quadratic penalty factor α-and then iteratively search for the optimal solution corresponding to the model, which can adaptively match the optimal center frequency and finite bandwidth of each intrinsic mode function (IMF) and realize the effective separation of the IMFs, and the general IMF expression is: Construction of the variational problem: Assuming that the original sequence f(t) is decomposed into K intrinsic IMFs, i.e., u k (t) with corresponding center frequencies and finite bandwidths and sparse, it is required that f (t) is equal to the sum of all u k (t) while pursuing the minimization of the sum of estimated bandwidths of each u k (t) i.e.: (1) where: k is the number of modes to be decomposed (positive integer), {u k } , {ω k } correspond to the kth modal component and the center frequency after decomposition, respectively, δ(t) is the Dirac function, * is the convolution operator, and e −jω k t is the phase-volume representation of the center frequency in the complex plane.
(2) Solution of the variational model.To find the optimal solution for the above variational model, the quadratic penalty factor α and the Lagrange multiplier λ are used to construct the augmented Lagrange expression, which transforms the constrained variational model into an unconstrained variational model, denoted as follows: where α can be used to ensure that the reconstruction accuracy of the signal places a limit on the bandwidth.Using the alternating direction multiplier method, {u k } , {ω k } and are updated alternately until the convergence condition is satisfied.
Stepwise decomposition of VMD Constructing correct and efficient training test samples is the key point to serve the real precipitation prediction.In other medium-and long-term prediction fields, the more commonly used sample construction methods that can serve the real world are 17 : Semi-stepwise decomposition (SSD), Fully stepwise decomposition (FSD), Single-model semi-stepwise decomposition (SSD), Single-model FSD, and Single-model FSD), Fully stepwise decomposition (FSD), Single-mode SSD (SMSSD), and Single-mode SSD (SMSSD).Single-mode SSD (SMSSD) and Single-model FSD (SMFSD).The first two methods require K (K is the number of modal components) models to be constructed simultaneously.The latter two require only a single model to obtain the final prediction results, providing higher operational speed.The sample construction processes for FSD and SMFSD are described below, respectively.
The important thing about the FSD technique is that the initial sequence of samples to be decomposed is only an initial sequence of the length of the number of response variables, and then the test period samples are subsequently appended one by one to the training set for decomposition.
(1) Sequence segments (S 1 , S 2 ,…, S m ) are decomposed to obtain K modal components, and the last m elements of each subsequence are extracted as explanatory variables.(2) The data S m+1 is appended to (S 1 , S 2 , …, S m ), and the extended sequence segments (S 1 , S 2, …, S m , S m+1 ) were decomposed, and the last m elements of each subsequence were extracted to continue to constitute the explanatory variables.By attaching new data one by one, the sequence segments are gradually expanded and the corresponding explanatory samples are extracted.For the last sample, the sequence segments (S 1 , S 2 , … S N-1 ), are decomposed and explanatory variables are extracted.(3) Response variables for all samples were extracted by decomposing the entire sequence segments (S 1 , S 2 , … S N ), obtained by delaying the values of the corresponding one time period, and finally all the samples were partitioned into training and test sets.
The biggest difference between SMFSD and FSD is that the precipitation results are directly used as the response variables instead of the results of each component, and the final summation and reconstruction are not needed, so the efficiency of the model operation is greatly improved.The construction process of SMFSD is shown in Fig. 1.
(1) The original precipitation series (S 1 , S 2 , … S N ) is divided into the training period (S 1 , S 2 , …, S p ) and test period (S p1 , S p+2 , …, S N ).(2) Sequence segments (S 1 , S 2 , …, S m ) were decomposed into K subsequences, and the last m elements of each subsequence were extracted as explanatory variables for the first sample, while the response variable was Sm + 1.Then, the sequence segments (S 1 , S 2 , …, S m , S m+1 ) is decomposed into K subsequences, and so on, until the sequence segment (S 1 , S 2 , …, S m , S m+1 , S N-1 ) are decomposed.Finally, N-m samples can be obtained.The first P-m samples are used as the training set and the remaining samples are used as the test set.Since the response variable is the actual precipitation data, SMFSD also requires only one model to be trained.

Modified AVOA
AVOA African Vultures Optimization Algorithm (AVOA) is an optimization algorithm proposed by Abdollahzadeh et al. 18 in, which has shown better optimization results in relevant tests compared to the Teaching and Learning Optimization Algorithm (TLBO), Gray Wolf Optimization Algorithm (GWO), Particle Swarm Algorithm (PSO) and Differential Evolutionary Algorithm (DE).The AVOA algorithm defines the objective as the vulture's hunger level: the worse the fitness value is, the vulture feels hungry and tries to search for food, which is the reality of the algorithm's optimization process.The AVOA iteratively proceeds as follows.
(3) (1) Phase 1.After randomly initializing the population position, the optimal or suboptimal individuals are selected according to the "roulette" rule for the next stage of the search.For the ith vulture in the population, the selection of its learning object can be expressed as follows: where, L is a user-defined parameter, located between (0,1), when L tends to 0, it contributes to the increase of population diversity, and vice versa, it accelerates the aggregation of the population; rand is a [0,1] uniformly distributed random number.(2) Phase 2: Define the starvation rate to realize the conversion of algorithm development and exploration process.Vultures start to search for food because of hunger, but if they are hungry, they don't have enough energy to fly for a long time over a long distance, and they can only search for food next to strong vultures to save energy.Therefore, the AVOA algorithm defines a starvation rate F as shown in Eq. ( 7), which increases as the iterative process progresses in order to facilitate the development process more likely.
where, h, z are uniformly distributed random numbers in the interval [−2,2], [−1,1] respectively; w is a user-defined parameter that controls the probability that the algorithm enters the exploration mode in the final stage; i iter is the current number of iterations in the course of the algorithm.(3) Phase 3: Spatial exploration.In this process, the AVOA algorithm uses a user-defined parameter P to determine what kind of exploration mode to enter, taking a value between (0,1).This phase is executed according to Eq. ( 8) to Eq. ( 9): where, P(i + 1) is the updated position of the vulture; X is [0,2] uniformly distributed random number.

Tent chaotic mapping initialization of populations
Excellent population initialization can speed up the convergence of the population and improve the performance of the system.The chaotic sequence generated by chaotic mapping has the characteristics of good traversal and strong randomness 19 , which can be used to replace [0,1] uniformly distributed random numbers to generate diverse populations.Therefore, choosing suitable chaotic sequences for population initialization is a key issue at this stage, and the use of chaotic sequences to initialize algorithmic populations is widely used at present 20,21 .
Tent chaotic mapping 22 with randomness, ergodicity, and regularity characteristics of chaotic sequences in large quantities, is one of the means commonly used to improve the representativeness of the initial population of optimization algorithms, and has been verified by examples to have better traversal than classical chaotic mapping 23,24 .The mathematical model of standard Tent mapping is as follows: Zhang et al. 25 considered that there might be small and unstable periodic points in the classical Tent series and introduced a random variable considering the population size N into Eq.( 17) to obtain the improved Tent mapping.
In using the above improved Tent mapping, it is found that the obtained chaotic variables are easy to go beyond the [0,1] boundary when the population size is small.Therefore, N = 100 is taken uniformly in this paper.
Chaotic mapping initialization population method According to the characteristics of Tent chaotic mapping, the steps of generating N chaotic individuals within the feasible domain are as follows.
Step 1 According to the number of dimensions D, take random initial values X 1i , X 2i , …, X Di , i = 1 in the interval [0,1].

Step 3
Stop iteration when i equals to the maximum number of iterations, each X sequence corresponds to each dimension to produce N chaotic individuals.

Selection of suboptimal vultures
In the standard AVOA algorithm, the first and second ranked vultures in the population in terms of fitness values will serve as the optimal and suboptimal vultures as potential learning objects for the population during the next iteration.However, the optimal and suboptimal vultures tend to be closer together in the later stages of the population iteration, resulting in the algorithm having been bound by a local optimum.Therefore, considering strengthening the ability of the AVOA to escape from local optimality at the later stage, this paper proposes (10) that the vulture with the third highest fitness ranking or even any vulture in the population can potentially be used as a learning object.On this basis there is the following mathematical model.rr= 0.9 + 0.1 * chaos(iter); Suboptimal vultures = second best individual for fitness value.else Strategy 1: Suboptimal vultures = any but the best and next best individual.end In the above suboptimal vulture selection model, chaos represents the chaotic sequence generated by the Tent mapping with the length of Maximum number of iterations; iter is the current number of iterations; and rand represents the random number generated by [0,1] uniform distribution.Whenever a suboptimal vulture is to be selected, a chaotic number rrgreater than 0.9 is first generated and then selected according to the model.Therefore, there is a lower probability of selecting vultures with weaker fitness values, and only vultures with the second highest fitness value are considered as suboptimal vultures for the initialization of the population using the Tent chaotic mapping.Thus, the exploration capability of the algorithm is improved without significant loss of exploitation capability.
In addition, the algorithm's own parameters also have a profound effect on the algorithm's performance.The standard African vulture optimization algorithm is parameterized with a step size of 0.1 for its four parameters L, P1, P2, and P3, resulting in a recommended parameter combination.In this study, a finer step size (0.05) was set for this set of parameter combinations, details of which can be found in Table 1.

DNC
DNC hybrid learning neural network, is a special kind of recurrent neural network with external memory (not trainable).At each time step t, a trainable controller decides the final output prediction by linearly combining the prediction information of the two parts after exchanging the information flow with the external memory based on the information flow at the moment t − 1.
Figure 2 illustrates the overall architecture of the DNC model.Input data is processed within the controller, and the generated memory information is written to the memory bank.The interaction between the controller and the memory bank is facilitated by the read and write heads. (

1) Controller η
The controller η in the DNC is a variant of the LSTM.at some time, the controller η obtains an input vector x t , Read R vectors r At each moment t, the controller η computes the output vector v t and an interaction vector ξ t , defined as follows: The controller η passes the information back by creating a loop in the computational graph, which in turn yields v t etc. 26 Finally, the η output vector is as follows.
(2) Read and write head and memory The controller η operates on the data in the memory through the read/write head.The read or write positions are determined by the corresponding weights, and the set of allowed weights at N positions is the non-negative quadrant of the standard simplex form in R N .
In the read operation, multiple read weights w r,1 t , w r,2 t , . . ., w r,R t ∈ � N are used to compute the weighted average of the contents, thus defining the read vector as follows.
In the write operation, the write weight w W t ∈ N is combined with the erase vector e t and the write vector v t to modify the memory storage matrix.
Here: ° denotes the Hadamard product; E is the All 1 matrix of N × M.
(3) Time-series memory link The temporal memory link holds information about the order in which the memory locations are written and is denoted as L t .L t i, j indicates the degree to which location i is written after location j is written, and each row and column of L t can be defined with location weights.
Here, p t is the priority weight and p t [i] denotes the degree to which position i was written last.
For read head i , define reverse weights f i t and forward weights b i t :

Model coupling
Step 1 K modal components are obtained by decomposing the precipitation sequence using the VMD algorithm according to different stepwise decomposition sample construction methods.
Step 2 In order to more accurately describe the number of pre-influence factors of each component, for the kth modal component according to the ACF and PACF to determine the corresponding influence lag month lag k taking Huaian as an example, 6 modal components are obtained, and the lag month of each component is lag k = [2,7,4,6,6,3].Step 3  Training and testing samples were generated using various stepwise decomposition methods, with a training-to-prediction sample ratio of 19:2, and normalization was applied based on the training samples.

Decomposition of data
The decomposition of the whole precipitation time series at one time will bring the information that can only be known in the future into the model training process, which will lead to the distortion of the prediction.However, due to the "boundary effect", decomposing the training set and test set separately often leads to poor prediction results, so the "stepwise decomposition" method is more scientific and reasonable.In order to better compare the prediction effects of "full decomposition" and "stepwise decomposition", the monthly precipitation series of the city of Huaian is used to analyze the decomposition effects of the two methods, and the analysis results are shown in Fig. 6.
In the actual precipitation prediction work, the past historical data belong to the training period, and the future results are the testing period.The first 227 samples of Huai'an precipitation series and the first 228 samples of the same station are selected, i.e., P (1:228) is obtained by adding one month of precipitation data to P (1:227) to simulate the actual precipitation prediction.
Using the VMD decomposition of P (1:228) and P (1:227) with the same parameters, subfigure (a) in Fig. 6 shows the first component in the respective decomposition results.P (1:228) and P (1:227) are identical except for the last sample, so the red and black lines in subfigure (a) should coincide, but they are separated at the allimportant right-hand boundary.Fundamentally, the red line represents the actual first component, i.e., the result of the full decomposition, while the black line represents the first component obtained in the actual application, which is clearly not equal.Therefore, from the point of view of additional sensitivity, the full decomposition method cannot reflect the actual application results, and subfigure (b) also reflects this difference.Extending this additional sensitivity to a longer test period, subfigure (c) further illustrates the decomposition difference between the full decomposition and the additional step-by-step decomposition, which can be seen to follow the entire test period.In summary, the "full decomposition" approach cannot serve the actual precipitation prediction.In addition, when the modal components are superimposed to obtain the reconstruction difference between the full decomposition and the stepwise decomposition in subfigure (d), it can be seen that the reconstruction result of the stepwise decomposition represented by the red line has a smaller difference with the actual precipitation and has a lower impact on the subsequent prediction.

Rainfall forecast
In order to further analyze the suitability of different decomposition methods for coupling with the prediction models, three sets of coupled models consisting of the "full decomposition", "FSD" and "SMFSD" methods of VMD are used to predict the monthly rainfall data of Zhengzhou city, Linyi city, Luan city and Huai'an city in 20 years.Next, the "full decomposition", "FSD" and "SMFSD" methods of VMD are used to predict the monthly rainfall data of Zhengzhou City, Linyi City, Liuan City and Huaian City in 20 years, so as to identify more scientific and reasonable coupled models with higher prediction accuracies.The rainfall from 2000 to 2018 was used as the training period, and the rainfall from 2019 to 2020 was used as the test period.The prediction results of different coupled models are shown in Fig. 7.
The figure shows the prediction results of the coupled models with different decompositions for the four cities, and it can be seen that the prediction accuracy of the three forecasting models is lower in the flood season of June-September, when the rainfall varies greatly.Among them, the rainfall in June-September 2020 in Liuan City changes drastically compared with the previous years, so the forecast results of the three coupled models here all have large errors.In contrast, the coupled model with SMFSD decomposition in subfigure (a) predicts the rainfall most closely with the actual rainfall; the coupled model with FSD decomposition in subfigure (b) can predict the overall trend of rainfall changes, but the forecasting accuracy is lower than that in subfigure (a) at the peak of rainfall, and the coupled model with full decomposition in subfigure (c) does not adequately explore and has a certain error in the cyclical trend of the rainfall time series during decomposition.The coupled model with full decomposition in subplot (c), due to the insufficient excavation of the cyclic trend of rainfall time series during the decomposition and the existence of certain errors, has a greater impact on the subsequent forecast, and can only roughly forecast the trend of rainfall, and the overall forecasting accuracy is obviously lower than that of the previous two.
In order to more intuitively show the prediction effect of the three different decomposition modes of the coupled model on each city, the linear fitting effect of the three is compared as in Fig. 8, and the average prediction accuracies of each model for the four city predictions are shown in Table 3.As can be seen from Fig. 9, the forecasting accuracy of the three coupled models is generally low for the city of Luan due to the presence of the sudden 2020 rainstorm in the city; compared to the other three cities, the forecasting effect is best for the city of Zhengzhou, which has a smoother trend in the rainfall time series.For all cities, the forecasting effect is SMFSD-MAVOA-DNC > FSD-MAVOA-DNC > VMD-MAVOA-DNC, Similarly, the results in Table 3 show that SMFSD-MAVOA-DNC has the highest average forecasting accuracy, and compared with the traditional VMD decomposition, the RMSE decreases by 7.42, the MAE decreases by 4.83, and the NSE increases by 0.05.NSE increased by 0.05, which fully indicates that the decomposition method of SMFSD is more suitable for actual rainfall forecasting.

Discussion
In order to further verify the prediction effect of the coupled models, four models, SMFSD-MAVOA-DNC, SMFSD-AVOA-DNC, SMFSD-MAVOA-LSTM, and MAVOA-DNC, are selected to conduct experiments in turn.Figure 9 shows the comparison of the prediction errors of each model for the four cities.
From Fig. 9, it is evident that the SMFSD-MAVOA-DNC model consistently demonstrates higher overall prediction accuracy under identical conditions.Additionally, it excels in predicting extreme values when compared to the other three models.Regarding the prediction performance for each city, it's observed that the models deliver poorer predictions for Liuan, a city with a larger standard deviation.Conversely, the prediction accuracy for Zhengzhou and Huaian, with smaller standard deviations, is relatively high.This indicates that the degree of data discretization in rainfall data has a noticeable impact on the prediction accuracy of each coupled model.However, the SMFSD-MAVOA-DNC model consistently maintains the highest performance, and the influence of data discretization on its prediction accuracy is minimal.
Table 4 demonstrates the results of rainfall prediction error comparison among the four cities, from which it can be observed that the SMFSD-MAVOA-DNC model performs the best in rainfall prediction in Zhengzhou, Linyi, Liuan and Huaian.Its average RMSE is 9.02, average MAE is 7.13, and average NSE is 0.94, which are much better than the other three models, showing higher prediction accuracy.In addition, the prediction errors of the four cities are relatively consistent, indicating that the model is more stable and representative for predicting rainfall in the Huaihe River Basin.
The combined model based on MAVOA-DNC is able to predict the general trend of the monthly rainfall series, however, compared with the SMFSD-MAVOA-DNC prediction model, its overall prediction effect is poorer, which fully demonstrates that the SMFSD decomposition effectively reduces the non-stationarity of the time series with high complexity and strong nonlinearity; meanwhile, the DNC combines the advantages of the recurrent neural network and computational processing that significantly improves the memory forgetting

Conclusion
Through the comparative study of different coupled models for rainfall prediction in Huaihe River Basin, we get the following conclusions: 1.For a single nonlinear and non-stationary rainfall time series, full VMD decomposition is not suitable for practical precipitation forecasting.After conducting a comparative study, it was determined that the best prediction accuracy can be achieved by using SMFSD decomposition.When compared to traditional VMD decomposition, the SMFSD coupled model demonstrated a significant improvement, with an average RMSE reduction of 7.42, MAE reduction of 4.83, and NSE increase of 0.05 in the four cities.This unequivocally demonstrates that the SMFSD decomposition method is better suited for real-world rainfall prediction.2. The systematic and valuable information obtained through SMFSD decomposition is applied to DNC for prediction, effectively leveraging DNC's strengths in prior learning of both long and short-term memory data.During the prediction process, multiple hyperparameters need to be configured, and the MAVOA algorithm introduced in this study provides superior optimization for machine learning models with a greater number of hyperparameters, thereby significantly reducing uncertainties in the modeling process.As a result, we achieved improved prediction results: the SMFSD-MAVOA-DNC model exhibited an average RMSE of 9.02, an average MAE of 7.13, and an average NSE of 0.94.When compared to the SMFSD-MAVOA-DNC and SMFSD-MAVOA-LSTM models, this model demonstrated the highest accuracy in rainfall prediction in the Huaihe River Basin. 3. In comparison with models proposed in other research articles, the SMFSD-MAVOA-DNC model adopts a more practical decomposition approach for rainfall prediction.It integrates multiple optimization and analytical techniques, enabling it to effectively address the challenges presented by nonlinear and non-stationary rainfall time series, while reducing uncertainties in the modeling process.This model provides a robust tool for rainfall forecasting.It is worth noting that, although the predictive performance of SMFSD-MAVOA-DNC is relatively good, there is room for improvement in predicting data with high levels of dispersion.Additionally, this paper analyzed the fitting capabilities of coupled models based on the VMD algorithm's full decomposition method.Other decomposition and denoising algorithms could be considered as subjects for future research, aiming to further enhance the accuracy and stability of rainfall prediction.

Figure 4 .
Figure 4. Location map of the study area.

Figure 8 .
Figure 8.Comparison of the linear fitting effects of the models.

Table 2 .
Monthly rainfall data summary table for various cities.

Table 3 .
Average prediction accuracy of each model.

Table 4 .
Comparison of relative errors of SMFSD-MAVOA-DNC with other models.